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ABSTRACT 


Rock fall is the act of releasing pieces of rock from bedrock slopes 
provoked by environmental factors, such as physical and chemical 
weathering, and vibration caused by rock blasting. The evaluation of 
the rate at which pieces of rock are released from the Quarry’s bench 
face through modeling in order to prevent possible accidents is 
necessary. The method used for this investigation involves taking 
field measurements of the slope height, slope angle, length and width 
of the bench. Rock samples were also taken from the field to the 
laboratory to carry out tests to determine the density, friction angle 
and slope roughness of the rock. With this, the models were made 
with rock fall program. The zone of the Quarry face prune to 
affectation by Rock fall was predicted to be 15 m from the bottom 
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crest of the IBD Quarry bench face that was occupied by rock pile 


and 6.5 m for a face without pile. Also, total kinetic, translational and 
rotational energies of the falling rocks were modeled at different 
points of the face; the bounce height and end-point were determined. 
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1. INTRODUCTION 

Rock fall is the act of releasing pieces of rock from 
bedrock slopes provoked by environmental factors, 
such as physical and chemical weathering, and 
vibration caused by rock blasting or heavy vehicular 
movement (Schumm and Chorley, 2017; Day et al., 
2018). Apart from the weathering rates, trigger 
mechanisms also determine whether Rock fall occurs. 
Rock fall can be activated in slope by freezing and 
thawing, geo-structural and seismic activities 
(Douglas, 2005, Gardner, 2017 and Bull et al., 2008). 
Wieczorek (2015) reported that human activities 
leading to decreased stability of hill slopes in hard 
rock and geological factors can also aid rock fall. 
(Selby, 2016). 


Rock fall can occur in different modes, such as 
freefall, where motion that occurs in the slope 
gradient is greater than 76° and the translation of the 
rock and the rotation of the block around its centre 
were observed (Ritchie, 2010; Azzoni et al., 2012) 
and air friction can minimally affects the velocity of a 
freefalling rock (Bozzolo and Pamini 2012). 
Moreover, at a reduced slope (gradient of about 45°) 
a bouncing rock gradually transform its movement to 
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rolling (Broilli, 2016). Hungr and Evans (2011) 
reported that a rolling rock is almost constantly in 
contact with the slope surface. Other type of motion 
over the surface of slope is sliding, which generally 
occurs at the beginning and ending stages of rock fall, 
the rock will usually stop because of frictional energy 
(Bozzolo and Panini, 2012). Therefore, the stopping 
of a falling rock depends on the mean slope gradient, 
the size of the rock and the material that covers the 
surface of the slope such as soil, scree and vegetation. 
Also, small rocks retard more easily than bigger 
rocks, because the total kinetic energy of small rocks 
is lower than that of bigger rocks during rock fall, 
trees can more easily stop small rocks and small rocks 
retard more easily in depressions between larger 
rocks on talus slopes (Kirkby and Statham, 2017). 


There exist empirical models, process-based models 
and GIS-based models for the Rock fall analysis. 
According to Keylock and Domaas (2010) empirical 
Rock fall models are referred to as statistical models 
and are based on relationships between topographical 
factors and the length of the run-out zone of one or 
more Rock fall events. 
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Tianchi (2013), carried out 76 major Rock fall test, 
established inverse logarithmic correlation 
relationship between the volume of the Rock fall and 
the ratio of the maximum vertical drop to the 
maximum horizontal distance travelled. Also the 
relationship of a positive logarithmic correlation 
between the volume of the Rock fall and the area 
covered by the fallen mass was formed. With this, 
Tianchi (2013) developed a model for an initial 
estimate of the extent of treating rock fall, if the 
volume can be estimated. Through comparative study, 
Moriwaki (2015) found a relationship between the 
angle of a line connecting the toe with the crown of 
the rupture, in the first instance, the ratio of the 
maximum vertical drop to the maximum horizontal 
distance travelled and, secondly, the landslide 
volume. Evans and Hungr (2011) on their side, 
suggest the Fahrbéschung principle to forecast run 
out zones of Rock fall events. The Fahrbéschung is 
the angle between a horizontal plane and a line from 
the top of a Rock fall source to the stopping point for 
any given rock fall. The line must follow the fall track 
of the boulder. On the other hand, Evans and Hungr 
(2011) suggested the minimum shadow angle as an 
alternative principle. It is the angle of a straight line 
between the peak point of the talus slope and the 
stopping point of the longest run out boulder. After 
investigating 16 talus slopes in British Columbia 
Evans et al., (2018) reported a minimum shadow 
angle of 27.5° which should only be used for a first 
approximation of the length of a Rock fall run out 
zone. Keylock and Dommas (2010) used Rock fall 
data to test three empirical models on their ability to 
predict the maximum length of Rock fall run out 
zones using simple topographic parameters. 


Process-based models describe the modes of motion 
of falling rocks over slope surfaces. Kirkby and 
Statham (2017) created a process based Rock fall 
model for transport of rocks that slide over slopes 
surface. This model calculated the velocity of the 
falling rock at the base of the cliff with the following 
expression: 


V S28 & DB vssccscccsices (equation 1) 


Where: 

v: velocity (ms~'); 

g: acceleration due to gravity (ms ~~) and 
h: fall height (m). 


Based on this velocity, the component of the fall 
velocity parallel to the slope surface was calculated, 
assuming that this component of the velocity is being 
conserved during the first impact of the rock on the 
slope surface. Then, the stopping position was 
calculated by the ratio of the fall velocity and a 
frictional force, which was determined by the 


dynamic angle of friction. Keylock and Domaas 
(2010) developed the simple dynamics Rock fall 
model, which is a process-based model based on the 
model of Kirkby and Statham (2017). The simple 
dynamics Rock fall model was tested using Rock fall 
data presented by Domaas (2008). Their model 
calculated the travel distance over the slope surface 
on the basis of the friction force according to Kirkby 
and Statham (2017) and the acceleration due to 
gravity. On the basis of calculated probabilities of 
modeled Rock fall travel distances, Keylock and 
Domaas (2010) concluded that the simple dynamics 
Rock fall model did not appear to hold a significant 
advantage over the empirical models tested in their 
study. In addition to the models of Kirkby and 
Statham (2017) and Key lock and Domaas (2010), 
there is a large group of process-based models that 
are rather similar (Wu, 2000; Bozzolo and Pamini, 
2012; Evans and Hungr, 2011; Bozzolo et al., 2015; 
Azzoni et al., 2012, Chau et al., 2017). Three factors 
correspond in all these models. First, these process- 
based models are two-dimensional slope-scale models 
that restricted falling boulders to move in a vertical 
plane. Consequently, lateral movements were not 
simulated. Secondly, the Rock fall track was defined 
as a composite of connected straight lines with a 
slope angle equal to the measured mean slope 
gradient on the represented segment of the Rock fall 
track as visualized. Finally, motions were simulated 
as a succession of flying phases and contact phases. 
The flying phase was simulated with a parabola 
equation based on the initial velocity in two directions 
and the acceleration due to gravity. The collision 
point of the rock on the slope surface was calculated 
with the intersection of the parabolic flying function 
and the straight slope segments. The first difference 
between these two-dimensional process-based models 
is that some of these models considered a falling rock 
with its mass concentrated in one point (Bozzolo and 
Pamini, 2012; Kobayashi and Cruden M.D 2018; 
Evans and Hungr, 2011; Azzoni et al., 2012), while 
other models considered bouncing, rolling and sliding 
as identical movements that can be described by a 
succession of impacts and bounces (Bozzolo et al., 
2015; Pfeiffer and Bowen, 2015). Models applying 
specific algorithms for calculating rolling and sliding 
velocities mainly used Coulomb’s law of friction, 


PF, = pp-m. g. cos@ equation 2 


Where, 

Fy is friction force (tangential to the slope surface) 
(kg/ms” 

uy is coefficient of friction; 

m is mass of the rock (kg); 

g is acceleration due to gravity (ms”); 
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© is mean slope gradient (°). 


The calculated friction force could then be used for 
calculating the sliding or rolling velocity of a rock 
after displacement over a given distance over the 
slope surface (Bozzolo and Pamini, 2012; Dijke and 
Westen, 2013; Kobayashi and Cruden M.D 2018; 
Evans and Hungr, 2011; Azzoni et al., 2012). Here, 
the friction coefficient is the most determining factor 
for the velocity. 


GIS-based models are those either running within a 
GIS environment or they are raster based models for 
which input data is provided by GIS analysis. GIS- 
based Rock fall models consist of three procedures. 
The first procedure identifies the Rock fall source up 
areas in the region of interest, the second determines 
the fall track and the third calculates the length of the 
runout zone (Hegg and Kienholz, 2016). Meissl 
(2014) developed two GIS-based Rock fall models 
using an empirical model for calculating the runout 
zone. The first model was Schattenwinkel. This 
model was based on the minimum shadow angle 
principle (Evans and Hungr, 2011). The second 
model of Meissl (2014) was called Geometrische 
Gefalle. This model was based on the angle of the 
shortest line between the top of the Rock fall source 
scar and the stopping point. Apart from these 
principles both models were identical, since both 
models used an identical module for calculating the 
fall track and the source areas (Callaghan and Mark, 
2013; Meissl, 2014). In a GIS-based distributed Rock 
fall model the terrain is represented by multiple 
rasters, which are derived from GIS data layers. Each 
raster represents a certain property of the terrain. 
Dijke and Westen (2013) developed a distributed 
model that performs a neighborhood analysis to a 


2. STUDY AREA 


DEM derived from an isocline map with an 
equidistance of 20 m, representing an area of 80 km’, 
A neighborhood analysis calculated the fall direction 
for each raster cell. The direction was determined 
towards the neighboring cell with the minimum 
height value. The velocity calculation of the model 
was based on an energy conservation principle as 
described by Scheidegger (2012). On the basis of this 
principle the velocity of the falling rock was 
calculated following, 


V =¥V, + 2g. (h —pyt 4x) equation 3 


Where, V is velocity of the falling rock (ms); 
Vois initial velocity of the falling rock (ms~'); 

g is acceleration due to gravity (ms); 

his fall height (m); 

uy is coefficient of friction; 

X is distance travelled over the slope surface (m). 


The velocity was calculated for each raster cell within 
the fall track, starting at a potential Rock fall source 
cell. Potential Rock fall source cells were defined by 
the mean slope gradient (>60°) and the value of the 
friction coefficient (u,) depended on the soil cover 
type. Dijke and Westen (2013) concluded that their 
model was able to compute the general distributions 
of Rock fall areas shown on their maps. The GIS- 
based distributed model of Meissl (2014) was called 
Sturzgeschwindigkeit and simulated two different 
modes of motion: freefalling and sliding. Meissl 
(2014) concluded that the simulated Rock fall 
patterns on high altitudinal zones generally 
correspond with the mapped patterns, but 
unfortunately no model accuracies were given. 


This study was carried out in IBD Impex limited Quarry, Afao, Ekiti state, located in longitude 5°19'32" E and 


latitude 7°42'25"N. 


\ > 
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Fig 1 Location Map of IBD Impex limited, Afao, Ekiti State, Nigeria (https://www.google.com) 
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3. MATERIALS AND METHODS 
The Quarry bench height was measured and the width was measured from the bottom crest of the bench under 
investigation to the top crest of the following bench. Also, slope angles for both quarries were determined by 
equation 4: 

Vertical distance 


Slope Angle = Tan™+ ————————— equation 4 
forizontal distance 


And the slope distance by equation 5 
Bench Hetght,H 


- »Where @ = angle of bench slope equation 5 
sin@ 


Slope distance, L = 
In order to estimate the weight of rock that will likely fall from the quarry face, density of the rock was 
determined in the laboratory by taking ten samples of intact rocks from the quarry site in a random manner, 
weight and of each sample was measured, then the measuring cylinder was filled with water to the desired level 
and the reading, Vo, was taken as initial volume. Separately, each of the rock samples was dipped inside of the 
cylinder and the reading Vj, as the final reading. The volume change was determined by the following relation: 


AV = V, — V, equation 6 


The average density of the sample was calculated by the expression in equation 7: 
= ightof s ri 
Average density = meh Tg SE Se ees equation 7 
Then the size of the likely boulders that can fall were estimated by measuring the length of the protruded edges 
of the piece of the rock, its volume was determined and consequently the weight of the possible falling rock was 
calculated by multiplying the density by the determined volume. Then the volume of the rock was multiplied 
with specific gravity to get the weight of the rock. 


Rock weight = V...SG, equation 8 
Where, V, is the volume of the rock 
SG, is the specific gravity of the rock 


3.1. MODELING 

The input parameters used modeling the rock fall at the face of the quarry are: angular velocity, coefficient of 
normal restitution, slope roughness, friction angle and coefficient of tangential restitution. In this study, the 
initial velocity is taken to be cero, because at this point the movement of the rock is very slow. On the other 
hand, the coefficient of normal restitution (Ry) and tangential restitution (Rr) were taken from the literature 
according to the case study of the rock that is similar to the rock under study, that is, Rn= 0.35 and Rr= 0.83. 


The Rock fall was modeled with the Rocfall program using the parameters measured which are slope height, 
slope angle, rollout distance and elevation of the rock. The slope was first created, the slope materials were 
defined and materials were assigned to the slope segments. The slope roughness was determined by tracing a 
straight a line of 1 m along the wavy face of the rock slope under study. The angle of inclination of the wavy line 
was measured at several points by sampling a normal distribution. The mean of the distribution which is equal to 
the original slope segment inclination and the standard deviation were calculated. 


4. RESULTS AND DISCUSSION 

Figure | shows the face of bench slope of IBD Implex Limited, revealed poor blasting operation, there appeared 
dangerously hanged intact rocks which could be unsafe to the loading and transportation crews due to possibility 
of rock fall. IBD quarry is largely made of granite rock. 
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Impex, Limited, Afao, Ekiti State, Nigeria. 


bench, IBD 
The data used to estimate the modeling of the rocks that were thrown are shown in Table | and Table 2. 
Table 1 Modeling data for IBD Impex, Limited 


2.79 | 2.64 |7.4] 12 | 15 | 90 37 0.35 0.04 0.83 0.04 1.5 4.6 
Table 2 Vertices of the IBD Impex limited quarry bench 
1 0 398 
2 15 398 
3 15 387 
4 30 38h 


The modeled face with blasted pile were shown in figure 2a and 2b with waving down sloping which is due to 
the zigzag positioning of the blasted rocks. When throw occurs, the seeder rock may likely hit any of the pointed 
area of the rock pile and this dynamism could change the direction of the seeder, cause the bounce height and 
end-point location to increase, even, the throw velocity could be made higher, thereby increasing the seeder 
throw energy. In this study, the bounce height of 0.5 m was recorded at the location 5 m from the bench lower 
crest, | m at location 10 m and from location 11.5 to 14 m, the bounce height ranges between 0.5 to 0.9 m. On 
the other hand, the rock end-point, where it finally came to a stop was 15 m from the bottom crest of the bench 
height (figure 2b). 


Bounce Height Envelope Horizontal Location of Rock End-points 


7 
\ 
! 


Height Above Slope (m} 
Number of Rocks 


+ 
10 20 
Location [m] Location [m] 


Fig. 2 The modeled IBD quarry bench face with the graph of bounce height and rock end point 
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TOTAL KINETIC ENERGY ENVELOPE Translational Kinetic Energy Envelope 


Fig. 3 The modeled IBD quarry face with the graph of the types of energies displayed 


Also, figures 3a, b, c, show different energies with which the rocks thrown were travelling: total kinetic energy 
was highest (900 J) at the initial point when seeder was released. As the rocks travelled the kinetic energy got 
reduced to 100 J at the location, 4 m from the bottom crest of the quarry bench, then it increased to 200 J at 
location 10 m and again descended to 60 J at 11 m, ascending to 100 J at 12 m and finally the kinetic energy got 
reduced to zero at 15 m; the rotational and translational energies follow the same waving pattern which range 
between 15 to 25 J within location 1 to 27 m, and between 50 to 500 J within 10 m from the lower crest of the 
bench respectively. The unstable manifestation of the energies was due to the waving nature of the rock mass 
pile at the quarry face. 


Bounce Height Envelope Horizontal Location of Rock End-points 


Fig.4. the rock throws bounce height and end-points at the face without rock pile. 


Moreover, figure 4a and b shows that the bounce height was at its maximum of 0.2 m at the 5 m location 
measuring from the quarry face and the end-point was at 6.5 m from the lower crest of the quarry face. 


Rotational Kinetic Energy Envelope Total Kinetic Energy Envelope 


Fig. 5 The modeled IBD quarry face without rock pile and its graph of energy types. 


In Mine design, safety of workers and equipment is paramount, especially during operations at the quarry face. 
In the light of this study, the end-point values will help to determine the distance that is safe for workers to stand 
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and the various types of energies determined will help to measure the impact of the Rock fall on the objects. 
When the quarry operations are to take place within lesser distance, protection must be provided against possible 
rock fall and the operation equipment must be constructed with integrity that is strong enough to arrest the 
destruction of the same by rock fall. In IBD Quarry, considering the height of Rock fall (10 m) analyzed, the 
end-point of 15 m means that the safe place for workers and equipment is any distance beyond 15 m from the 
bottom crest of the Quarry when rock pile is still on the face site and beyond 6.5 m in the case where the face is 
free of rock pile. Also, with the knowledge of the energy of the rock fall, the likely impact can be estimated and 
the design of the equipment could be such that the effect of the impact would be reduced to bearable minimum. 


5. CONCLUSIONS 

In this study, the zone of the Quarry face prune to 
affectation by Rock fall was predicted to be 15 m 
from the bottom crest of the IBD Quarry bench face 
with rock pile and 6.5 m for a face without pile. Also, 
total kinetic, translational and rotational energies of 
were modeled at the different points of the face; the 
bounce height and end-point were determined. 


5.1. RECOMMENDATION 

This study is recommended for the Open Pit Mine 
project design department of Mining Industries. But 
caution should be made to carry out all necessary 
process of data collection and laboratory test 
according to Industrial standard before application. 
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